High dimensional behavior of the Kardar-Parisi-Zhang growth dynamics 
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We investigate analytically the large dimensional behavior 
of the Kardar-Parisi-Zhang (KPZ) dynamics of surface growth 
using a recently proposed non-perturbative renormalization 
for self-afHne surface dynamics. Within this framework, we 
show that the roughness exponent a decays not faster than 
a ~ 1/d for large d. This implies the absence of a finite upper 
critical dimension. 
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The study of the non equilibrium dynamics of rough 
surfaces and interfaces has received a great deal of at- 
tention in the last years HQ. Both theoretically and 
experimentally many efforts have been devoted to single 
out the traits and features shared by apparently different 
phenomena. In this context, by analogy with equilibrium 
statistical mechanics, the search for universality classes 
is a central task. The Kardar-Parisi-Zhang || equation 
(KPZ) is, for surface growth, the main contribution in 
this direction. It is the minimal Langevin equation cap- 
turing the essence of many different growth models be- 
yond the Gaussian linear theory HI!- It reads 



dh(x, t) 
dt 



v^ 2 h+ ^(V/i) 2 + T)[x,t). (1) 



where h(x, t) is the surface profile, x is the position in a 
d-dimensional substrate, r\ is a Gaussian white noise, v 
and A are constants. The KPZ equation also describes 
the behavior of directed polymers in random media |Q, 
systems with multiplicative noise J5|, and it is related to 
the Burgers equation ||. 

A central quantity of interest is the roughness W(L) 
of a system of linear size L, defined as 



W 2 {L) = ±- d Y,[h{x,t)-hf 



(2) 



where h = {1/L )^2 x h(x,t). In many seemingly unre- 
lated growth processes the large scale properties of the 
roughness are observed to be scale invariant and univer- 
sal; i.e. in the stationary state W(L) ~ L a and corre- 
lations decay on a typical time t s ~ L z , with universal 
exponents a and z. These critical exponents are not in- 
dependent, as a consequence of the Galilean invariance 
of the related Burgers equation fj]^J a + z = 2. It is thus 
sufficient to focus the attention on one exponent, say a. 



The theoretical analysis of the KPZ is extremely diffi- 
cult. Apart from the d — 1 case, where a special symme- 
try makes an exact solution possible with a — 1/2, the 
situation is still quite controversial despite the large effort 
devoted to the problem. In particular, the fundamental 
issue of the existence of an upper critical dimension d c , 
above which the exponents recover their mean-field (or 
infinite dimensional) values (a — 0) @], is highly debated 
H . The application of field theoretical tools presents an 
inherent problem: one indeed finds that the fixed point 
controlling the rough phase of the KPZ is not accessible 
to perturbation expansion in A; this fact renders stan- 
dard field theoretical tools inadequate for this problem. 
Early applications of non-perturbative methods such as 
functional renormalization group |^] and Flory-type ar- 
guments |l(J suggested that d c = 4, in agreement with 
a l/d-expansion ]ll[ | around the d = oo limit. Later 
the mode-coupling approximation led to contradictory re- 
sults suggesting the existence of a finite d c jl2| or d c = oo 
Jl3| . Arguments for a finite d c based on directed |l4| or 
invasion jnj percolation have also been proposed. More 
recently a detailed analysis of a d = 2 + e perturbative 
expansion revealed a singularity at d — 4 p6| , leading 
Lassig to the conclusion that d c — 4 is the upper critical 
dimension of the KPZ dynamics |Tij ]. 

Numerical simulations of models in the KPZ universal- 
ity class markedly disagree with this last conclusion ||, 
showing that a > at least up to d = 7 (l8| . In particu- 
lar numerical results suggest a large-d behavior a ~ 1/d 
in agreement with early conjectures J|, p. 75]. Both of 
these conclusions were confirmed by a recently proposed 
renormalization group (RG) approach Jl9| . The key idea 
of this approach is that the geometric scaling of the grow- 
ing surface can be ascribed to a scale invariant dynamic 
process, which builds the same correlations at all length- 
scales. This scale invariant dynamics is the fixed point of 
the RG transformation, which is derived by consistency 
requirements of the description of the same system at 
two different scales. Analogous ideas, implemented via a 
real space RG, have proved to be quite powerful to inves- 
tigate the critical properties of non-equilibrium, strong 
coupling problems |2(|. The implicit nature of the RG 
transformation, which is similar in spirit to the idea of 
phenomenological RG [|TJ, allows us to avoid the use of 
hierarchical lattices, a source of incontrollable approxi- 
mations, specially in high dimensions. Remarkably, the 
exponents predicted by the RG are in excellent agreement 
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with numerical simulations up to d = 7. 

In this Letter we analyze the large-d behavior of this 
RG approach and show that it predicts that the rough- 
ness exponent a vanishes not faster than l/d for d > 1. 
This rules out the existence of a finite upper critical di- 
mension. In what follows we expose the essential con- 
cepts of this method and apply it to the analytical study 
of the KPZ dynamics in the large-d limit. 

Consider a growing surface, whose dynamics, at the 
microscopic scale, is defined in terms of a stochastic equa- 
tion, such as Eq. (1), or by a discrete model. If we par- 
tition the d + 1-dimensional space in cells of lateral size 
Lfe = 2 k Lo and vertical size h k , we obtain a static descrip- 
tion of the surface at the coarse grained scale L k : With 
some majority rule each block is declared to be empty or 
filled. For each substrate cell i the number h(i) of filled 
blocks on top of it identifies the interface configuration, in 
units hk, at scale Lk- Note that h k is an independent pa- 
rameter of the static description. Scale invariance implies 
that if hk is properly chosen, the coarse grained system 
looks similar at all (large enough) length-scales L k . The 
optimal geometric description, which best exhibits scale 
invariance, in our case, is that with hk oc W(Lk) ~ L% of 
the same order of typical height fluctuations over a dis- 
tance Lk- In the RG procedure, we shall fix hk = 2W(Lk) 
in order to have a scale invariant description of the sur- 
face (see ref. [H for details). The coarse-graining pro- 
cedure, which defines the static description in terms of 
blocks of size Lk, also induces a flow of the microscopic 
dynamics towards an effective dynamics at the same scale 
Lk] this is defined in terms of the transition rates for the 
addition of an occupied block. The main feature of KPZ 
dynamics is lateral growth J3pj|, and this suggests the 
following minimal parametrization of the growth rates at 
the generic scale Lk is 

r[h(i) -> h(i) + l] = l+.T fc ^ max[0, h(j) - h{i)]. (3) 

jnni 

The first term is the contribution of the vertical growth 
(i.e. random deposition) and the second term is the con- 
tribution of lateral growth. Indeed the sum over neighbor 
block sites j counts the area of the vertical surface ex- 
posed towards site i. Xk is then the ratio of lateral to 
vertical growth rates. We shall come back later, in the 
conclusions, to the approximations implied by Eq. (||). 

In order to derive the RG transformation, let us par- 
tition a system of size L into £ d cells of size Ljl. We 
observe that the roughness W 2 (L) can be written as the 
average roughness W 2 (L/£) inside single cells plus the 
fluctuations of the average height among different cells. 
The second contribution, within the description at scale 
Ljl, is simply given by the roughness uj 2 (£, x) of a system 
of I cells - with x being the dynamic parameter at scale 
L/t - times the square height of the cells [2W(L/£)] 2 . 
Hence we find 



W 2 {L) = W 2 (L/£) [1 + Aoj 2 (£, x)] 



(4) 



which is the basis of the RG approach. With I = 4, 
L = L k+2 and x = x k , it gives W 2 (L k+2 ) /W 2 {L k ). 
The same quantity can be alternatively computed using 
Eq. (^) with 1 = 2, once with L = Lk+2 and x = x k +i, 
and a second time with L = L k +\ and x = x k . The 
consistency of the two calculations yields an implicit RG 
transformation 



l + 4w 2 (2,.x fe+1 ) 



1 + 4w 2 (4, Xk) 
l+4co 2 (2,Xk) 



(5) 



for the dynamic parameter x k ■ The attractive fixed point 
x* = lhrifc^oo Xk (if it exists) identifies the scale invariant 
dynamics and Eq. (Q) with Xk = x* finally yields the 
roughness exponent 



a = lim log 2 



'W 2 (L k+ i) log[l + 4w 2 (2,x*)] 



W 2 (L k ) 



2 log 2 



(6) 



Eqs. (^|, |6|) are the starting point of our analysis. A more 
detailed discussion of their derivation can be found in 
Rcf . jlj| . We note here that the existence of an attractive 
fixed point x* implies that the process is "self-organized" : 
No fine tuning is necessary in order to observe the critical 
behavior. 

A key observation is that, since lo 2 (£,x) — > for x — > 
oo, x* = oc is a fixed point of Eq. (j^) corresponding to 
a = 0. Therefore the RG scheme allows, in principle, 
for the occurrence of a finite upper critical dimension d c 
(a = for d > d c ) and the existence of a finite attractive 
fixed point for all d is a non-trivial prediction. A finite 
stable fixed point was found in Ref. for d = 1, . . . , 8 
using Monte Carlo methods to compute uj 2 (£,x). The 
same method was also applied to the Gaussian theory 
[A = in Eq. (1)], recovering the result d c = 2, i.e. a = 
for d > 2 0j. Though very powerful, the Monte Carlo 
method cannot be pushed to very high dimensions nor 
does it provide an explicit analytic behavior of a as a 
function of d. 

In the following we study analytically the large-<i limit 
of the RG in order to extract its predictions on the ex- 
istence of a finite upper critical dimension and on the 
large-d behavior of the roughness exponent. 

The technical difficulty lies in the explicit calculation of 
the functions uj 2 (£,x) for I = 2, 4. For d> 1 we expect 
a <C 1, which means that surface fluctuations u>(£, x) ~ 
£ a ~ 1 + a In £ + . . . are of order 1 . This suggests that for 
a system of small size £ we can reasonably account for 
the fluctuations of the interface if we allow h(i) to take 
only two values: h(i) = ho or h(i) = ho + 1. This drastic 
approximation has the advantages of making the explicit 
computation feasible on the one hand, and of providing 
a lower bound for the exponent a on the other. We shall 
come back later to this important issue. Let us only 
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stress, for the time being, that a lower bound on a is 
sufficient to exclude the existence of a finite d c . 

In the above approximation, growth can only occur on 
"low" sites (h(i) = ho). This means that Eq. (^) is only 
valid if h(i) = ho and the rates vanish on "high" sites 
(h(i) = ho + 1). It is convenient to classify the possible 
configurations {h(i}} by the number n of "high" sites. 
The roughness Eq. (Q) of each configuration of n "high" 
sites is the same and is equal to (1 — nj£ d )nj£ d and the 
dynamics involves only transitions from configurations 
with n to configurations with n + 1 "high" sites. We can 
then group all configurations{/i(i)} with n "high" sites in 
the same effective state with a great simplification of the 
structure of the master equation (the state with n = £ d 
is equivalent to the flat surface n = 0). The only non- 
vanishing transition rates r{n — ► n+1) are obtained from 
Eq. (Q) summing on all possible final configurations and 
taking the average on the initial configurations, which 
leads to 



2 d+1 ln2. 



r(n 



I) 



n + xfl r 



(7) 



The first term here accounts for vertical growth, which 
can occur only on the £ d — n "low" sites. The second 
term is the contribution of lateral growth and is the 
average number of lateral walls (i.e. the surface between 
"low" and "high" sites) in configurations with n "high" 
sites. Assuming that "low" and "high" sites are randomly 
distributed, each "low" site has on average 2dn/£ d "high" 
neighbor sites and therefore 



a n ~ 2d(£ d 



(8) 



The distribution of "high" sites is actually not random 
but we have verified numerically that, for large enough 
dimensions, Eq. (^) provides a reasonable approximation 
p3| . Combining Eqs. (0, ||) one easily obtains the prob- 
ability p n of state n in the stationary state of the master 
equation 



Pn = Po 



r(0 



1) 



1, 



(9) 



where po 

Z^n=0 Pn ' 



P0 = 1 



r(n —> n + 1) 

is fixed by the normalization condition 
1. A simple calculation leads to 

i 

1 + 2dxk 



2dxk 



2d In ^ + In 



£ d + 2dx k 



(10) 



The roughness of configurations with n particles, using 
Eq. (||), is (1 — n/£ d )n/£ d which, averaged over the dis- 
tribution p n (as specified by Eq. @ and (10)) gives 

u 2 %x)=po^— (11) 
zdxk 

where we have assumed dxt 3> 1 and £ d ^> 1. Combining 
Eq. ( |Tl| ) with the RG equation (JsJ) we obtain, to leading 
order in d, a fixed point 



(12) 



consistent with the assumption dxk 3> 1. Using now 
Eq. (JsJ) it is straightforward to find, to leading order in 
d > 1, 



1 1 



3(ln2) 2 d' 



(13) 



Furthermore, we can also analyze the stability of the fixed 
point. The derivative of the RG transformation Xk+i = 
R(xk) of Eq. (H), at the fixed point, is 



R'(x* 



-1 



1 1 



2 In 2 d 



o(d~ 2 ) 



(14) 



Since |i?'(x*)| < 1 we can conclude that the fixed point 
is attractive Vd Therefore we find a finite, stable fixed 
point x* with an exponent a > for all d, which is the 
main result of this Letter. This excludes the occurrence 
of a finite upper critical dimension, d c , which would show 
up, in the present framework, in a stable fixed point at 
x* = oo for d > d c . 

Let us now discuss the validity of the approximations 
used. We neglected configurations with h(i) > ho + 2 or 
equivalently deposition processes on a "high" site. The 
rate of this process, on a state with n "high" sites, is 
^up(^) = n. Our approximation is then valid if r(n — > 
n + 1) 3> r up (n). This condition fails when the process 
is close to complete a new layer, i.e. for n ~ £ d . More 
precisely the deposition on "high" sites is not important 
for 



n > 1 > 



2dx' 



(15) 



Since x* ~ 2 d , the approximation is correct for £ — 2 
V<£ Fig. [l] shows the that the approximation to the 
LHS of Eq. (|) 1 + Auj 2 {£ = 2,x) is good already for 
d = 7. The approximation is much less accurate for £ — 4 
and, as a consequence, fluctuations in the system of size 
£ = 4 arc underestimated. This means that our approach 
underestimates the RHS of Eq. (||) and consequently also 
its value at the intersection point with the LHS. This 
value is directly related to the roughness exponent by 
Eq. (|(|) and therefore the restriction of height fluctuations 
leads to a lower bound to the exponent a. Accordingly - 
since the LHS of Eq. (||) decreases with x - Eq. (jlj) gives 
an upper bound to the true fixed point parameter x* . 
Fig. |l| illustrates this analysis for d — 7. Fig. || shows a 
comparison of the present analytical estimates [Eqs. jl^ ) 
and (||)] and the results of ref. g, jj]§ and |TJ. 

Besides the approximations of the present calculation, 
which, as we have argued, provide a lower bound to a, 
it is also worth discussing the approximations of the RG 
method itself. In this respect we observe that Eq. (^) is 
a minimal parametrization of the dynamics, in the sense 
that it allows for the minimal proliferation in the RG cap- 
turing the relevant features of KPZ growth. In principle, 
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more proliferation parameters can be included in order 
to improve the accuracy of the method. It is important 
to note, however, that the range of typical fluctuations 
h(i) — h(j) ~ l a allowed in the RG calculation is small 
and the one-parameter approximation in Eq. (|3|) to the 
scale invariant dynamics is reasonable. This is confirmed 
by the accuracy of the RG predictions in finite dimensions 
p9[ and it is expected to improve as a — > 0. Therefore the 
inclusion of additional proliferation parameters in Eq. (|^) 
is not expected to change the nature of the fixed point 
and of our main conclusions. Let us also point out that 
usually small cells analysis becomes very accurate in high 
dimensions. An extension of the RG procedure to cells 
of larger size, going beyond the present approximations, 
provides in principle a systematic way to improve our 
prediction which is currently under investigation |22f . 

In conclusion we have shown that the recently pro- 
posed []l9| real space RG predicts that the roughness 
exponent a decreases not faster than 1/d as d — > oo 
[Eq. Jl3|)]. This implies that there is no finite upper 
critical dimension in the KPZ universality class and it 
suggests that theoretical arguments leading to d c = 4 
should be reconsidered. 

We acknowledge interesting discussions with A. Mari- 
tan, G. Parisi, A. Stella, C. Tebaldi and A. Vespignani. 
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"High" sites rather than being randomly distributed, 
tend to coalesce in droplets. This means that the ap- 
proximation is worst for n ~ l d /2. This effect is not very 
important for small sizes I in d S> 1. 
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FIG. 1. Graphic analysis of Eq. (5) from the present ap- 
proximation and from Monte Carlo evaluation for d — 7. 




FIG. 2. Value of a from the present calculation (full line) 
and from simulations of ref. 8 and 18 for d — 1, ... ,7. In view 
of the approximations involved (see text) we obtain a lower 
bound. Inset: fixed point value x*(d) vs d. The theoretical 
prediction (full line) is an upper bound to the true x*(d). 
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